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^ Abstract 
O 

Continuum reduction and Monte Carlo simulation are used to calculate the heavy quark potential 
and the string tension in large TV" Yang-Mills theory in four dimensions. The potential is calculated 
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out to a separation of nine lattice units on a lattice with extent six in each direction. 
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I. INTRODUCTION 



Through dimensional transmutation, the dimensionless classical coupling constant of 
QCD is transformed to a running coupling controlled by a physical scale Aqcd- Roughly 
speaking, weak coupling perturbative approximations to processes with momentum transfer 
Q are expansions in l/ln(Q/ Aqcd)- This gives excellent results when Q/Aqcd is large. 
When one cannot rely on that, an alternative is the large N approach. The gauge group is 
generalized to SU(N), and the expansion is in 1/N. 

Although the large N limit of Yang-Mills theory is still out of reach, it is known that it 
enjoys several simplifications. Among them is the possibility of reducing the space-time vol- 
ume without affecting certain physical results pp. Continuum reduction [2l|3j, i.e. reduction 
to a physical size of order 1/ Aqcd avoids some of the difficulties [U El E] with reduction to a 
single space-time point. The addition of double trace terms to the action [7] is an alternative 
approach. 

It is now practical to obtain good results by combining continuum reduction with numer- 
ical simulations needing only modest resources. In previous numerical work [HI El QUI ED] , 
it has been shown that not only bulk quantities but also physical quantities based on Wil- 
son loops are accessible. The previous results based on Wilson loops have been in three 
dimensions. In the work described here, we show that the method is still practical in four 
dimensions. 

We have used continuum reduction and Monte Carlo simulation to calculate the heavy 
quark potential and the string tension in large N Yang-Mills theory in four dimensions. An 
important aspect of the method is that reduction allows the calculation of infinite volume, 
infinite N Wilson loops that are larger than the reduced lattice. In particular, in this work, 
the heavy quark potential is calculated out to a separation of nine lattice units on a lattice 
with extent six in each direction. The results for the string tension are compatible with 
those obtained on large lattices at smaller N. 

II. METHODS 

The standard Wilson Yang-Mills action with gauge group SU(N) is used. In the large 
N limit, g 2 is taken to zero with the inverse 't Hooft coupling b = -h^ held fixed. 
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We report results with N = 37, 47, and 59 on a lattice of size 6 4 . This is large enough so 
that in the range of available b, the lattice is not too coarse but is still small enough to give 
a manageable computational cost. 

In four dimensions, the useful range of couplings and physical lattice spacings a on a 
given lattice is more limited than in three dimensions. For 6 4 , the system becomes unstable 
to the bulk transition for small b. The smallest b we have used is 0.3450, which is just 
above the unstable point. For sufficiently large b, the center symmetry breaks in one lattice 
direction, and reduction no longer holds in that direction. This is the large N limit of the 
finite temperature phase transition [21 H2] ■ For 6 4 it occurs at about b = 0.3515. Our 
calculations are at b = 0.3450, 0.3480, and 0.3500 with most of the results at 0.3480. As 
measured by the critical size for the finite-temperature transition, that is a range from about 
L c = 4.4 to 5.6. Thus the possibilities for testing scaling are quite limited. For this lattice 
size, this entire useful range of b is in the region that is metastable to the bulk transition. 
Nevertheless, we encountered no bulk transitions during the simulations. 

As described in previous work [3J, we update the gauge field configurations with heat 
bath and over-relaxation methods. In this work, one update will mean one heat bath sweep 
followed by one over-relaxation sweep. Measurements of the Wilson loops were separated 
by ten updates. The values for the Wilson loops are based on 1200 measurements. 

The measurements of Wilson loops are made on smeared configurations. The use of 
smeared links improves the measurement of Wilson loops. They enhance the overlap of the 
space-like sides of the Wilson loops with the ground state. This increases the signal relative 
to the fluctuations and simplifies the t behavior of the loops [TB]. The smearing is a four- 
dimensional version of the method used in [10]. One lattice direction is arbitrarily chosen 
as the "time" direction. Links in the remaining three spatial directions (but not in the time 
direction) are smeared. After the Wilson loops in these "time" -space planes are measured, 
the process is repeated with each of the other lattice directions chosen as "time." 

When smearing the links in spatial directions, only staples in spatial planes are used. 
One step in the iteration takes one from a set Uj?\xi,X2,t) to a set Uj? +1 \xi, X2, t). Before 
reunitarization, the weight of uj:\xi, X2,t) is (1 — /) while that of each staple is //4. There 
are two parameters, namely, the smearing factor / and the number of smearing steps n. We 
use / = 0.45 and n = 5 so that t — fn — 2.25 and the associated length scale is y/r = 1.5. 
The use of finer smearing steps / = 0.1, n = 25 or a larger length scale / = 0.45, n = 10 
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was more costly and did not lead to further improvement. 

Data were collected on planar, rectangular Wilson loops of size k x j with k and j ranging 
from 1 to 9 and with j the extent in the "time" direction. The j decay of the loops is fit to a 
simple exponential. To avoid a possible distortion from a combination of smearing and very 
small error bars at the shortest separation, loops that are 1 x j and fcxl are not included 
in the fit. The rate of the exponential decay is taken as the static quark potential at the 
separation k. The k dependence of the exponential is fit to obtain the string tension. With 
all quantities in lattice units, the three parameter potential that is used in the fits is 

m{k) = aa 2 k + c + -j- (1) 

We expect Co to be positive and of order 1/6 and c\ to be negative and of order 1/6 or 0(1). 

Errors in all quantities at a fixed b and N are obtained by jack knife with single elimina- 
tion. 



III. RESULTS 



We have results for N = 47 at b values of 0.3450, 0.3480, and 0.3500. In addition, there 
are results for iV = 37 and JV = 59 at b = 0.3480. 

An example of fits to Wilson loop data is given in Fig. [TJ which is for b = 0.3480 and 
N = 47. At a fixed k, the decay in j of loop data is fit to the form Ae~ m ^ k ^ using k and 
j from 2 through 9, inclusive. This gives the potential m(k) at separation k which is then 
plotted as a function of k in Fig. [2] A fit of the potential for the b = 0.3480 and iV = 47 
data to the form of Eq. [T] gives 

aa 2 = 0.099 ± 0.016 (2) 

To verify that N is sufficiently large, we have the results for N = 37, 47, and 59 in Fig. [3} 
The string tensions from fits to the N = 47 and N = 59 data agree. 

For a check of scaling in the limited range available, we can compare the potentials for 
b = 0.3450, 0.3480, and 0.3500 all with N = 47. (For b = 0.3450, it was necessary to restrict 
the largest dimension of the loops to 8 to obtain useful fits.) The relative physical scales 
are set from ratios of the critical lengths L c (b) at which the center symmetry breaks in one 
direction. These are determined from the results in [3]. With 

LJb) 



we have s = 4.4/5.2, 1, and 5.6/5.2 for b = 0.3450, 0.3480, and 0.3500, respectively. In 
Fig. |4j we plot ms verses k/s and see that there is rough agreement of the physically scaled 
potentials. 

IV. CONCLUSION 

A comparison with results on large lattices is in order. This can be done by comparing ei- 
ther bare quantities at the same tadpole improved coupling or by comparing a dimensionless 
ratio. Our result for the string tension in lattice units is aa 2 = 0.099 ±0.016. This is at bare 
coupling b = 0.3480 which corresponds to a tadpole improved bj = &(□) = 0.182. (The aver- 
age plaquette (□) is normalized to approach one in the weak coupling limit.) As it happens, 
Lucini, Teper, and Wenger [HI [TS] have a large lattice, N = 8 result aa 2 = 0.116 ± 0.001 at 
the same &/. This falls within our much larger uncertainty range. 

An appropriate dimensionless ratio of physical quantities is T c j ^fa where, in our case, 
aT c = 1/L C . At bj = 0.182, our result is T c j ^fo = 0.61 ± 0.05. After extrapolations to 
infinite volume, infinite N, and zero lattice spacing, Lucini, Teper, and Wenger [151 EE E] 
obtain 0.597 ± 0.004. 

In conclusion, we have used large N continuum reduction and Monte Carlo simulation 
to calculate the heavy quark potential and the string tension in Yang-Mills theory in four 
dimensions. The results are compatible with those obtained on large lattices at smaller N. 
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FIG. 2: The static quark potential m(k). 
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FIG. 3: The N dependence of the potential. 
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FIG. 4: Scaling behavior of the potential. 



10 



